MZ-TH/11-41 
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We present next-to-leading order corrections in the leading colour approximation for jet rates in 
electron-positron annihilation up to seven jets. The results for the two-, three-, and four-jet rates 
agree with known results. The NLO jet rates have been known previously only up to five jets. The 
results for the six- and seven-jet rate are new. The results are obtained by a new and efficient 
method based on subtraction and numerical integration. 
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INTRODUCTION 

Precise calculations for multi-parton final states are 
extremely important for an accurate description of the 
physics at the LHC experiments. The precision and accu- 
racy is reached by including next-to-leading order (NLO) 
QCD corrections in the theoretical prediction. However, 
multi-leg NLO calculations are far from trivial. In par- 
ticular the virtual corrections for processes with many 
external legs have been considered to be a challenge for 
a long time. In the past years there has been tremen- 
dous progress both in the refinement of the traditional 
Feynman diagram approach and in the development of 
the unitarity method. With the currently available tech- 
niques one is able to obtain results for processes involv- 
ing six- or seven-point one-loop functions. Examples of 
such calculations are the NLO corrections to the pro- 
cesses pp -> V + 3, 4 jets [1-3], pp -> WW + 2 jets [4], 
pp -> tt + 2 jets [5], pp ->• ttbb [6, 7] or e + e~ -> 5 jets 
[8]. Most of these calculations are based on the unitar- 
ity method. The results involving one-loop seven-point 
functions have been obtained only in the leading-colour 
approximation. 

In this letter we report on results obtained by a com- 
pletely different method. We calculate the NLO correc- 
tions by numerical methods [9-15]. To this aim we ex- 
tend the subtraction method from the real emission part 
to the virtual part. This renders the integration over 
the loop momentum finite and we perform the loop inte- 
gration together with the phase space integration in one 
Monte Carlo integration. In this letter we demonstrate 
the potential of this method. We calculate the NLO cor- 
rections for the process e + e _ — > n jets in the leading 
colour approximation up to n = 7. 

The LEP experiments have collected data on the jet 
rates in electron-positron annihilation. This data to- 
gether with a theoretical calculation can be used to ex- 
tract the strong coupling a s . The precise value of a s is 
an important input parameter for all LHC experiments. 

The results for the two-, three-, and four-jet rates have 
been known for a long time and we use them to verify 
the validity of our method. The NLO corrections to the 
five-jet rate have been calculated recently. The results on 



the six- and seven-jet rate are new. Note that the NLO 
corrections to the seven-jet rate require the calculation 
of a one-loop eight-point function. To the best of our 
knowledge this is the first time that physical observables 
depending on a one-loop eight-point function have been 
calculated. 

We have written a numerical programme based on re- 
cursion relations. This programme can in principle cal- 
culate the NLO corrections to n-jet observables in the 
leading colour approximation for all n. Changing the 
number n amounts to changing a single line in the nu- 
merical programme. Within our method we have a very 
moderate growth in the CPU time per evaluation as n 
increases. The practical limitations arise from the fact 
that the number of evaluations required to reach a cer- 
tain accuracy increases with n. This behaviour is already 
present at the Born level and not inherent to our method. 

The application of our method is not limited to 
electron-positron annihilation. The calculation of pro- 
cesses in hadron-hadron collisions requires only minor 
modifications and will be addressed in the near future. 



THE NUMERICAL METHOD 

In this section we briefly outline the method for the 
numerical computation of NLO corrections. More details 
can be found in [9]. 

In electron-positron annihilation the contributions at 
leading and next-to-leading order for an infrared-safe n- 
jet observable O are given as 



(0) L0 = J O n da 
(O) 



(1) 



NLO = / O n+1 da* + J O n da v . 

n+l n 

Here a rather condensed notation is used. da B denotes 
the Born contribution, whose matrix elements are given 
by the square of the Born amplitudes with (n + 2) par- 
ticles. Similar, da K denotes the real emission contribu- 
tion, whose matrix elements are given by the square of 
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the Born amplitudes with (n + 3) particles, da v gives 
the virtual contribution, whose matrix elements are given 
by the interference term of the one-loop amplitude with 
(n + 2) particles, with the corresponding Born amplitude. 
The individual contributions at next-to-leading order are 
divergent and only their sum is finite. We rewrite the 
NLO contribution as 

(0) NLO = (0)NLp + (0) NLO ai + (0) NLO io ^ (2) 

In order to render the real emission contribution finite, 
such that the phase space integration can be performed 
by Monte Carlo methods, one subtracts a suitably chosen 
piece [16-19]: 
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This is the standard application of the subtraction 
method. Within our method we use the same idea for 
the virtual part: 



(O)virtual = 2/ #„ Re J 
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Here, G^ le denotes the integrand of the bare one-loop 
amplitude: 
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^soft' ^coii an< ^ ^uv are l° ca l subtraction terms, approx- 
imating the soft, collinear and ultraviolet singularities of 
the integrand. The subtraction terms have to be added 
back, which is done with the term 
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(O) 



NLO 
insertion 



O n (I + L) ® da 1 



(6) 



The insertion operator I contains the integrated subtrac- 
tion terms from the real emission part, the insertion op- 
erator L contains the integrated subtraction terms from 
the virtual part plus the ultraviolet counterterm from 
renormalisation. By construction, all three contributions 
(O)^p, (0)NLp al and <0>^° tion are finite. It should 
be noted that all quantities arc defined at the amplitude 
level and are calculated using recurrence relations. This 
allows for a fast implementation. 

Let us discuss eq. (4) in more detail. The subtraction 
terms ensure that the result of the loop integration is fi- 
nite. This docs not yet imply that we can perform the 
loop integration numerically by integrating the four loop 
momentum components from minus infinity to plus infin- 
ity. The integrand still has singularities for certain values 
of the loop momentum, i.e. if one or more of the propa- 
gators in the loop goes on-shell. This can be avoided by 



deforming the loop integration into the complex space, 
which is the second essential ingredient of our method. 
We devised and implemented an algorithm for the con- 
tour deformation in two variants. One variant is based on 
an auxiliary Feynman parametrisation [9, 14], the other 
one works directly in loop momentum space [12]. The 
two variants yield identical results, with the direct loop 
momentum deformation method being more efficient for 
higher jet rates. 



JET RATES IN ELECTRON-POSITRON 
ANNIHILATION 

In this paper we consider jet rates in electron-positron 
annihilation. The jets are defined by the Durham jet 
algorithm [20]. For the Durham jet algorithm the reso- 
lution variable is defined by 



Vij = 



2mm(E?, E)) (1 - cos %) 

Q 2 • 



(7) 



where Ei and Ej denote the energies of particles i and j, 
and Oij is the angle between the three-momenta pi and 
Pj. Q is the centre-of-mass energy. 

The production rate for n-jet events in electron- 
positron annihilation is given as the ratio of the cross 
section for n-jet events divided by the total hadronic cross 
section 



0"n-jct(«) 



(8) 



crtot(w) 

The arbitrary renormalisation scale is denoted by u. The 
jet rates can be calculated within perturbation theory. 
The perturbative expansion of the jet rate reads 

Rn(ji) = (9) 



2tt 



2tt 



In practise the numerical programme computes the quan- 
tities 



Cn-jct(u) 

<to(aO 
2tt 



2ir 



(10) 



B n (M) + 0«), 



normalised to cro, which is the LO cross section for 
e + e~ — > hadrons, instead of the normalisation to cr to t- 
From the expansion of the total hadronic cross section 

Ctot 

a {n)(l + \c F ^-+O{al)\ (11) 



0-tot(j«) 



we obtain the relations between the coefficients A n , B n 
and the coefficients A n , B n : 
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B n — B r , 



:C F A n . 



(12) 
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It is sufficient to calculate the coefficients A n and B n for a 
fixed renormalisation scale which can be taken conve- 
niently to be equal to the centre-of-mass energy /i = Q- 
The scale variation can be restored from the renormali- 
sation group equation 
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The values of the coefficients A n and B n at a scale /i are 
then obtained from the ones at the scale fio by 

A„(/x) = A n (/j|)), 

S n (Ai) = B„( A/0 ) + fc^/3 ln^^A„( A/0 ).(14) 

We can expand the perturbative coefficients A n and 
B n in 1/JV C : 
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In this article we calculate the coefficients A nt \ c and B n ^\ c 
for n < 7. 



NUMERICAL RESULTS 

In this section we present the numerical results. We 
calculate the leading order coefficient A rh \ c and the next- 
to-leading order coefficient B n ^\ c for n < 7 for the fixed 
renormalisation scale equal to the centre-of-mass energy 
Ho = Q. We take the centre-of mass energy to be equal 
to the mass of the Z-boson Q = mz- The calculation is 
done with five massless quark flavours. 

We first compare our approach with well-known results 
for two, three and four jets [21, 22]. Figure (1) shows a 
comparison for the two-jet rate for the NLO coefficient 
S2,ic between the numerical calculation and an analytic 
calculation. It should be noted that for the two-jet case 
the dependence on y cut enters only through the real emis- 
sion contribution. We observe an excellent agreement. 
The corresponding comparisons are shown for the three- 
jet rate in figure (2) and for the four-jet rate in figure (3). 
Again we observe an excellent agreement. The results 
for the five-, six- and seven-jet rates are given in table (I) 
for a few selected values of the jet parameter y cut . Note 
that a reasonably large value for the higher jet rates re- 
quires a small jet parameter y cut . 

Let us have a closer look at the performance of our 
method. Figure (4) shows the CPU time required for one 
evaluation of the Born contribution, the insertion term 
in eq. (6) and the virtual term in eq. (4) as a function 
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FIG. 1: Comparison of the NLO corrections to the two-jet rate 
between the numerical calculation and an analytic calculation. 
The error bars from the Monte Carlo integration are shown 
and are almost invisible. 
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FIG. 2: Comparison of the NLO corrections to the three- 
jet rate between the numerical calculation and an analytic 
calculation. The error bars from the Monte Carlo integration 
are shown and are almost invisible. 



of the number of final state partons n. One clearly sees 
that the insertion term is almost as cheap as the Born 
contribution. But more important is the fact that for 
all contributions the CPU time per evaluation increases 
only very moderately as a function of n. Within our 
method all three contributions should scale asymptoti- 
cally as n A [23]. Fitting a polynomial function n a to the 
numbers shows that this is indeed the case. In principle 
we could further improve the n 4 scaling behaviour to v? 
by replacing the four-gluon vertex by two effective three- 
valent vertices [24]. We emphasise that the virtual part 
has the same scaling behaviour as the Born contribution. 
This moderate growth imposes almost no restrictions on 
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FIG. 3: Comparison of the NLO corrections to the four-jet 
rate between the numerical calculation and an analytic calcu- 
lation. The error bars from the Monte Carlo integration are 
shown and are almost invisible. 



2/cut 






0.002 
0.001 
0.0006 


(5.0529 ± 0.0004) • 10 3 
(1.3291 ±0.0001) • 10 4 
(2.4764 ± 0.0002) ■ 10 4 


(4.275 ± 0.006) ■ 10 5 
(1.050 ±0.026) • 10 6 
(1.84 ±0.15) • 10 6 


2/cut 


*c A 

—Ae,ic 


N o R 


0.001 
0.0006 


(1.1470 ±0.0002) • 10 5 
(2.874 ±0.002) • 10 s 


(1.46 ± 0.04) • 10 7 
(3.88 ±0.18) • 10 7 


ycut 




<R 


0.0006 


(2.49 ± 0.08) ■ 10 6 


(5.4 ±0.3) ■ 10 8 



TABLE I: Perturbative coefficients for the five-jet rate, six-jet 
rate and seven-jet rate. 



the number of final state partons to which our method 
can be applied. The practical limitations arise from the 
fact that the number of evaluations required to reach a 
certain accuracy increases with n. This behaviour is al- 
ready present at the Born level and not inherent to our 
method. Altogether, the calculation of the seven-jet rate 
takes a few days on a cluster with 200 cores. 

CONCLUSIONS 

In this letter we reported on NLO corrections in the 
leading colour approximation for jet rates in electron- 



FIG. 4: CPU time required for one evaluation of the Born 
contribution, the insertion term and the virtual term, respec- 
tively, as a function of the number of the final state partons 
n. The times are taken on a single core of a standard PC. 



positron annihilation up to seven jets. The calculation 
is based on a new and powerful method, where the loop 
integration is performed numerically. To the best of our 
knowledge this is the first time that physical observables 
depending on a one-loop eight-point function have been 
calculated. We are planning to extend this method to 
hadron-hadron collisions and to include all terms in the 
colour expansion. 
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